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ABSTRACT 

I present a review of incompressible magnetohydrodynamic (MHD) turbu- 
lence in a strongly magnetised plasma. The approach is phenomenological even 
where a more rigorous theory is available, so that a reader armed with paper, 
pencil and some determination may be able to work through most of the physics. 
The focus is on the inertial-range spectra for very large (fluid and magnetic) 
Reynolds numbers. These theories of the inertial-range are built on two impor- 
tant facts: (i) Kraichnan's insight that the turbulent cascades are a result of 
nonlinear interactions between oppositely directed wavepackets of Elsasser fields; 
(ii) these oppositely directed wavepackets do not exchange energy, but contribute 
only to changing each other's spatial structures. I begin with a description and 
critique of the Iroshnikov-Kraichnan theory, and explore the fundamental depar- 
tures necessitated by the anisotropic nature of the turbulence. Derivations of 
the inertial-range spectra of four regimes of MHD turbulence — the balanced 
weak, balanced strong, imbalanced weak and the imbalanced strong cascades 
- are then presented. The need for studying the spectra of imbalanced turbu- 
lence when the waves on the outer scale have a short correlation time is briefly 
discussed. 

Subject headings: MHD — turbulence 



1. Introduction 

Astrophysical systems offer numerous settings for turbulence in magnetised plasmas, 
such as the Sun, the solar wind, accretion discs, the ionised interstellar medium, molec- 
ular clouds and clusters of galaxies (Armstrong, Rickett & Spangler 1995, Horbury 1999, 
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Biskamp 2003, Kulsrud 2005, Vogt & EnBlin 2005, Schekochihin & Cowley 2007 and McKee 
& Ostriker 2007). Magnetohydrodynamics (MHD) offers the simplest description of the dy- 
namics of a magnetised plasma. MHD is a single fluid description of a conducting magnetised 
plasma. Sources for the magnetic field, B(x,t), could be external or due to currents in the 
plasma. Lorentz forces act on the fluid and affect the evolution of the velocity field, v(x,t). 
If fluid motion is subsonic, incompressibility is a good approximation and the system can be 
described by the two vector fields, B and v. 

Stirring a magnetised plasma on spatial scale L creates velocity and magnetic field 
perturbations whose nonlinear evolution gives rise to a turbulent cascade to smaller spatial 
scales and is ultimately dissipated by viscosity and resistivity into heat. Naturally occuring 
processes can be somewhat different from stirring; for instance, the energy sources could be 
distributed in some way in space. In the case of the solar wind, the ultimate source of energy 
for MHD waves is the Sun, whereas in the interstellar medium the energy sources could be 
supernova remnants located in the spiral arms of our Galaxy. In these cases, MHD waves 
are generated in some regions of space, from which they propagate away, interact and make 
the plasma turbulent. Modelling this inhomogeneous situation is not easy. So, to simplify 
the physics of a difficult problem, we make the standard assumption that we are dealing 
with a statistically homogeneous and stationary situation. Thus stirring is included as a 
source term in the Navier-Stokes equation for the velocity field, with appropriate statistical 
properties. Moreover, stirring is indispensible in numerical simulations, which complement 
analytical studies. 

This review aims to be a self-contained description of MHD turbulence. However, there 
could be a couple of places where the reader may have to refer to other sources, and I 
hope I have made this clear in the text. This review could also seem somewhat one-sided 
in that much of it is based on work by me and my collaborator^. There are, of course, 
other points of view and the interested reader should consult the literature; these could be 
Iroshnikov (1963), Kraichnan (1965), Shebalin, Matthaeus & Montgomery (1983), Goldreich 
& Sridhar (1995, 1997), Ng & Bhattacharjee (1996), Cho & Vishniac (2000), Biskamp & 
Miiller (2000), Maron & Goldreich (2001), Cho, Lazarian & Vishniac (2002), Galtier et 
al. (2000, 2002), Lithwick & Goldreich (2003), Miiller, Biskamp & Grappin (2003), Galtier, 
Pouquet & Mangeney (2005), Boldyrev (2005), Miiller & Grappin (2005), Lithwick, Goldreich 
& Sridhar (2007), Beresnyak & Lazarian (2008), Chandran (2008), Perez & Boldyrev (2009) 
and Podesta & Bhattacharjee (2009). 



1 A topic on which I have not worked is imbalanced weak MHD turbulence, and the account in § 6 is 
based on my reading of Lithwick & Goldreich (2003). 
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In this review the focus is on the inertial-range spectra of the magnetic and velocity 
fields in a strongly magnetized incompressible fluid, in the limit of very large (fluid and 
magnetic) Reynolds numbers. § 2 begins with the equations of incompressible MHD which 
govern the dynamics of the magnetic field, B(x,t), and the velocity field, v(x,t). Both 
stirring and dissipation are included in the equations but will be dropped later. This is 
because in many cases they do not directly influence the inertial-range spectra for very large 
Reynolds number^. Alfven and Slow waves emerge as linear perturbations of a uniformly 
magnetised plasma, and their polarisations and dispersion relations are discussed. Nonlinear 
waves are introduced and their role in enabling a broad physical picture of MHD turbulence 
is sketched. § 3 discusses the Iroshnikov-Kraichnan theory of MHD turbulence and the 
problems affecting this theory. Attempts at remedying the defects lead to new physics and 
the next few sections are devoted to the theories that have emerged over the last fifteen years. 
However, I have dispensed with the chronological in favour of a more logical presentation: 
thus § 4 and 5 are on balanced weak and strong cascades, respectively, followed by imbalanced 
weak and cascades in § 6 and 7, respectively. In § 8 I briefly discuss the need for studying the 
spectra of imbalanced turbulence when the waves on the outer scale have a short correlation 
time. § 9 offers comments on matters that have not been discussed in the text. 

2. The physical system and MHD turbulence 
2.1. Incompressible MHD 

The equations describing incompressible MHD are given by, 



V ■ B = 0, V -v = (1) 

d t B - V x x B) = r]V 2 B (2) 

( B V") B 

d t v + (y ■ V) v = - Vp + - — + vV 2 v + f (3) 

47rp 

where p is the uniform and constant mass density of the fluid, r\ is the resistivity, and v 
is the kinematic viscosity. f(x,t) is the rate of velocity stirring, assumed incompressible, 
V- / = 0. In many cases of interest to problems in turbulence, / is stochastic, and its 



2 Dissipation does play an important role in the case of imbalanced weak MHD turbulence and this is 
discussed in § 6. Stirring can affect the spectra in the imbalanced strong cascade and this is discussed in § 8. 
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statistics is assumed to be given, p is the ratio of the total pressure to the mass density, 
which is determined by requiring V* v = 0. Thus we need to solve a Poisson equation, 



V 2 p = V- 



B V)B ■ 

- (v-V)v 



(4) 

with appropriate boundary conditions to determine p. Note that, in this incompressible 
limit, the thermodynamic orgins of p are not so apparent. Equations (0Q) — (j4j), together with 
initial and boundary conditions, completely determine the evolution of the vector fields, 
B(x,t) and v(x,t). Equivalently, the system can also be described by the Els as ser fields, 
W ± , defined by 

W + = v JL, W- = v + -J= (5) 

When p is uniform, we have V • W ± = 0. It turns out to be very useful to describe 
incompressible MHD turbulence in terms of the Elsasser fields, instead of the magnetic and 
velocity fields. 

Below we ignore explicit consideration of stirring and dissipative terms, because our 
focus is on the description of the inert ial-range at high Reynolds numbers. 



2.2. Transverse waves 

Uniformly magnetised static solution: 

v = 0, B = B z (6) 

In the uniformly magnetised state the Elsasser fields, = ^Vaz, where Va = B / ^/Anp 
is the Alfven speed. When this state is perturbed, we can write, 

W + = -V A z + w + , W = V A z + w . (7) 

It must be noted that the perturbations are not necessarily small. The perturbations obey 
the equations 



(d t + V A d z ) w + + (w ■ V) w + = -Vp 
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(d t - V A d z ) w- + (w+ ■ V) w~ = -Vp (8) 
where p is determined by requiring that V- w = 0. The linearised equations are 

(d t + V A d z )w+ = -Vp 

(d t -V A d z )w- = -Vp (9) 

If we consider solutions that vanish at infinity, V 2 p = implies that p = constant for the 
linearised problem, and we have linear waves with 

w + = F+(x,y,z-V A t) 

w = F-{x,y,z + V A t) (10) 

where F ± are arbitrary vector functions with V- F ± = 0. It is clear that w + is a wave that 
translates in the positive z-direction with speed V A , and w~ is a wave that translates in the 
negative ^-direction with speed V A ; wave propagation is along the magnetic field lines. We 
can restate this in fc-space by assigning frequencies, c t ; ± (fc) to the waves with dispersion 
relations, 



cj + (k) = V A k z , u~{k) = -V A k z (11) 

Being transverse waves, we must have k- F ± (k) = 0, where F^ is the Fourier transform of 
F^ . In other words, for any given k, the waves have two independent degrees of polarisation, 
constrained to lie in the plane perpendicular to k. The component perpendicular to z is called 
the Alfven wave, and the other orthogonal component is called the Slow (magnetosonic) wave. 



2.3. Nonlinear solutions and the general nature of MHD turbulence 

It is remarkable fact that equations (|S]) admit exact, nonlinear solutions: it may be 
verified that, when either w + or w~ is equal to at the initial time, the nonlinear terms, 
(w~ ■ V) w + and (w + ■ V) w~ , vanish for all time. Therefore there are exact, nonlinear 
solutions of the form, 
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F ± 



(x, y,zT V A t) , 



w 



T _ 







(12) 



In other words, wavepackets that travel in one direction (i.e. either in the positive or negative 
z-direction) do not interact nonlinearly among themselves. Nonlinear interactions occur 
only between oppositely directed wave-packets. From this fact, Kraichnan (1965) came to 
an important conclusion: 

(i) Incompressible MHD turbulence in a strongly magnetised plasma can be described as being 
the result of (nonlinear) interactions between oppositely directed wavepackets. 

It is worth noting another important result that we can prove by manipulating equa- 
tions (jSD, namely 



This result is true even in the presence of nonlinearity, and is equivalent to the conserva- 
tion of total energy and cross-helicity. Our second important conclusion now follows from 
equation ( {TBI : 

(ii) Wavepacket collisions conserve ±ve energies separately, and collisions can only redis- 
tribute the respective energies in k-space. This is true regardless of how imbalanced the 
situation is. For instance, if stirring puts much more energy into +ve waves than —ve 
waves, the +ve waves cannot transfer any part of their energy to —ve waves. They can only 
scatter the —ve waves and redistribute them (and vice-versa) in k-space. 



It is very useful to begin with an account of a theory of incompressible MHD turbulence 
in a strongly magnetised plasma due to Iroshnikov (1963) and Kraichnan (1965), even if 
current views hold it to be incorrect. The physical picture is more transparent in the latter 
work and our description below follows this. Assume statistically steady, isotropic excitation 
of ±ve waves, both with root-mean-squared (rms) amplitudes wl <C V a , correlated on 
spatial scale L, which can be referred to as the stirring scale or outer scale. Nonlinear 
interactions between oppositely directed wavepackets create ±ve wavepackets on smaller 
spatial scales. So let us consider the nature of collisions between wavepackets of scale A C L, 
where A is still much larger than the dissipation scale. We will assume that the most effective 
collisions occur between wavepackets of similar sizes, an assumption that is also referred to 




(13) 



3. Iroshnikov— Kraichnan (IK) Theory 
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as locality of interactions. From equations ([8]) it follows that, in one collision between a +ve 
packet and a — ve packet, the rms amplitude of either packet is perturbed by an amount, 



Sw x ~ -f (14) 

We imagine that a +ve wavepacket of size A goes on to suffer a number of collisions with 
— ve wavepackets of similar scale, and vice versa. Successive perturbations add with ran- 
dom phases, so that in n collisions the rms amplitude of the perturbation will increase to 
n l / 2 {w\/VA)- Therefore, the number of collisions for perturbations to build up to order unity 
is 



w 2 



N >~\t) (15) 

The cascade time is the time taken for collisions to build upto order unity. Since each 
collision takes time ~ X/Va, the cascade time is, 



tx ~ N x — (16) 

We now use Kolmogorov's hypothesis that the flux of energy cascading through the inertial- 
range is independent of the particular spatial scale belonging to the inertial-range. i.e. we 
assume that the energy flux, 



7* (17) 



is independent of A. This implies that in the inertial-range, 



/ A \ 1/4 

w x ~ w L ( — J (18) 

Thus we have determined the A-dependence of the rms amplitude in the inertial-range. The 
same result can also be expressed in fc-space in terms of the (isotropic) energy spectrum, 
E(k)k 3 ~ w\, where k ~ 1/A: 



(19) 
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However, we are not quite finished with our task. In our estimates, we assumed that 
N\ ^> 1 and we need to verify if this is true. Using equations (|15p and (|18p . we can estimate 
that N\ oc \- 1/2 . Therefore N x increases as A decreases; in other words, N\ increases as 
we go deeper into the inertial-range. So, it looks as if our phenomenological theory of the 
IK cascade is self-consistent, and we might be tempted to conclude that the inertial-range 
spectrum of MHD turbulence is given by equation (jT9l) . However, the physics is not so 
simple. 

3.1. Problems with the IK theory 

When N\ 3> 1, nonlinear interactions between oppositely directed wavepackets are weak, 
and we should be able to calculate better than we have done above. There is a well-developed 
theory, called weak turbulence, that describes weakly nonlinear interactions between nearly 
linear waves — see Zakharov, L'vov & Falkovich (1992). In this approach, perturbation 
theory is used to derive kinetic equations in fc-space describing wave-wave interactions. 
Then the inertial-range energy spectra emerge as stationary solutions carrying a flux of 
energy to large k. The lowest order of perturation theory, corresponding to what are called 
3-wave interactions in weak turbulence, is what is needed to make rigorous the IK cascade. 
Here, the rate of change of energy (or wave action) at a given k is due to either (i) the 
coalescence of two other waves with wavevectors ki and k 2 , or (i) decay of k into two waves 
with wavevectors k\ and k 2 . For any given k, all k\ and k 2 contribute, so long as they satisfy 
certain resonance conditions that depend on the dispersion relations satisfied by the linear 
waves. In the case of incompressible MHD, these 3-wave resonance conditions are somewhat 
special, because of the existence of the exact, nonlinear solutions given in equation (fT2"j) . Since 
there are nonlinear interactions only between oppositely directed wavepackets, it means that 
ki and k 2 must belong to oppositely directed wavepackets. Then the 3-wave resonance 
conditions may be written as, 

kx + k 2 = k 3 , uf + uj^ = uj^r or (20) 

where ujf = iV^A^, uj 2 = TVAk 2z , lo^ = ±VAk z and = TVAk z . Note the resemblence 
to "momentum-energy" relations in quasi-particle interactions in quantum field theory or 
condensed matter physics. These resonance conditions have a remarkable property, noted 
by Shebalin, Matthaeus & Montgomery (1983), namely that one of k\ z or k 2z must be zero. 
This fact may be readily seen by considering the ^-component of the wavevector resonance 
condition, k\ z + k 2z = k z , together with the frequency resonance condition. If, for instance, 
we choose ki and k to be + waves and k 2 to be a — wave, then the frequency resonance 
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condition is, k\ z — k^ z — k z . In this case, it is clear that = 0, implying that k\ z = k z . 

Hence waves with values ofk z not present initially cannot be created by nonlinear interactions 
during wavepacket collisions. In other words, there is no cascade of energy to smaller spatial 
scales in the z -direction, and hence the turbulence must be anisotropic. 

The somewhat surprising fact is that, attempting to create an anisotropic version of the 
IK theory, leads us very far from the IK theory and to entirely new physics. Before we get 
into details, it is useful to introduce some terminology. When the random stirring on the 
large scale L generates +ve and — ve waves with equal rms amplitudes (or energy input), 
we will refer to the situation as balanced; otherwise, the turbulence is called imbalanced. 
The case of imbalanced turbulence is also refererred to as turbulence with non zero cross 
helicity, or non zero (v - b), where b is the magnetic field perturbation. Also, the turbulence 
will be called weak when perturbation theory is applicable (when the cascade time is larger 
than the collision time); otherwise the turbulence will be called strong. Thus there are four 
cases to consider: balanced & weak, balanced & strong, imbalanced & weak and imbalanced & 
strong. Weak turbulence, whether balanced or imbalanced, can be described rigorously using 
the kinetic equations of weak turbulence. However strong turbulence, balanced or imbal- 
anced, does not permit such a rigorous description and we are limited to phenomenological 
descriptions and numerical simulations. 



4. Balanced Weak MHD Turbulence 

The cascade is due to the resonant 3-wave interactions discussed above; see Ng & 
Bhattacharjee (1996), Goldreich & Sridhar (1997), Galtier et al. (2000) and Lithwick & Gol- 
dreich (2003). Consider again steady, balanced and isotropic excitation of ±ve waves, with 
amplitudes t«i < 14 and scale L, which initiates a weak cascade to smaller spatial scales. 
As we have discussed above, there is no transfer of energy to smaller spatial scales in the 
direction parallel to the mean magnetic field. No parallel cascade implies that wavepackets 
of transverse scale A -C L have parallel scales L. 

It turns out that the different polarisations of the Alfven waves and Slow waves play 
a very important role in the turbulent cascade. Viewed physically, we may think of the 
(w^ 1 ■ V) terms as being responsible for the nonlinear cascadqj. During collisions be- 



3 This statement is not entirely accurate, because the Vp term also contributes a quantity of the same 
order of magnitude which does not have this precise form. However, the argument advanced in the text 
suffices for our phenomenological account in this review. 
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tween two oppositely directed wavepackets, we may think of (to 1 * 1 • V) as coming from the 
wavepacket that offers a perturbation, and the w T wavepacket as the one that suffers the 
perturbation. For A <C L, the polarisation of an Alfven wave is such that its w ± is very 
nearly perpendicular to z, whereas w ± in a Slow wave is almost parallel to z. Then we 
estimate that (w ± ■ V) ~ (w\/\) when the perturbing wavepacket is an Alfven wave, and 
(w^ 1 ■ V) ~ (w\/L) when the perturbing wavepacket is a Slow wave. Clearly an Alfven wave 
is a much stronger perturber than a Slow wave. Thus, it comes as no surprise that Slow 
waves are mostly scattered by Alfven waves, but do not significantly perturb Alfven waves; 
the cascade of Slow waves is slaved to the Alfven waves. Hence in our estimates below we 
ignore Slow waves, and restrict attention to Alfven waves. 

A collision between a +ve Alfven wavepacket and a — ve Alfven wavepacket now takes 
time ~ L/Va, because the wavepackets have parallel scales ~ L. In one collision a wavepacket 
is perturbed by amount, 



The perturbations add with random phases. The number of collisions for perturbations to 
build up to order unity is 



5w x 



(21) 




(22) 



The cascade time is 




(23) 



Kolmogorov's 



hypothesis of the A-independence of the energy flux 




(24) 



implies that in the inertial-range, the rms amplitudes are given by, 




1/2 



(25) 
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The energy spectrum, E(k), is highly anisotropic in fc-space, and now depends on both k z 
and k±. The precise dependence of E on k z is not so interesting, because there is no parallel 
cascade; it suffices to know that E is largely confined to the region \k z \ < Lr 1 . However, 
the dependence of E on kj_ is of great interest, because this has been established by the 
cascade in the transverse directions, described above. Accounting for this anisotropy, we can 
estimate, E(kj_, k z ) k\ L^ 1 ~ w\, where k± ~ A -1 . Now, using equation fl25|) . we have 

E(k ± ,k z ) ~ for \k z \ <L- 1 (26) 

Being a theory based on weak turbulence, the turbulent cascade can be described rigor- 
ously using kinetic equations for energy transfer; see Galtier et al. (2000) and Lithwick & 
Goldreich (2003). 

Having obtained the inertial-range energy spectrum of the balanced weak cascade, we 
need to perform checks of self consistency regarding the assumed weakness of the cascade. 
Using equations ff22|) and (T25]) . we estimate that 



V A \ 2 A 



iV A ~ — j (27) 

We must first verify that the cascade initiated at the stirring scale L is weak to begin with, 
which implies that we must have Nt 3> 1- From equation (|27|) . we see that this can be 
satisfied if wl <C Va- In other words isotropic, balanced stirring initiates a weak cascade 
if the rms amplitudes at the stirring scale is much less than the Alfven speed. A more 
geometric way of stating this is that a weak cascade is initiated if stirring bends field lines 
(of the mean magnetic field) only by small angles. 

Equation (j27p implies more interesting consequences: the transverse cascade strengthens as A 
decreases, because N\ decreases when A decreases. Therefore the assumption of the weakness 
of interactions must break down when N\ ~ 1, which happens when A ~ A* ~ L(wl/Va) 2 - 
Therefore, the inertial-range spectrum of equation (T26]) is valid only for transverse scales 
larger than A*. 



5. Balanced Strong MHD Turbulence 

When JVa ~ 1, the assumption of the weakness of the cascade is no longer even approx- 
imately valid. The turbulent cascade turns strong, and was first described by Goldreich & 
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Sridhar (1995). They conjectured that N\ remains of order unity for smaller values of A all 
the way down to the dissipation scale. This is equivalent to assuming that the cascade time 
remains of order the wave period, a condition that may be referred to as critical balance. 
The balanced weak cascade described above turns into a balanced strong cascade for A < A*. 
Instead of following the transition from weak to strong, we prefer to describe the inertial- 
range of the balanced strong cascade when it is initiated at the stirring scale itself because 
there is less clutter in the description. 

Consider steady, balanced and isotropic excitation of ±ve waves, with amplitudes »l ~ 
Va and scale L, which initiates a strong cascade to smaller spatial scales. Wavepackets of 
transverse scale A < L can now possess parallel scales that are smaller than L. This is 
because the resonance conditions are not in force when nonlinear interactions are strong so, 
in addition to a transverse cascade, a parallel cascade can also occur. Let eddies of transverse 
scale A possess parallel scale Aa, which is as yet an unknown function of the transvarse scale 
A. Critical balance implies that 



A A A 

t x ~ — ~ — (28) 
w x V A 

As before we invoke Kolmogorov's hypothesis of the A-independence of the energy flux 



— (29) 



which implies that in the inertial-range 



w x ~V A [^y\ A\ ~ L 1 ^ 3 A 2 ^ 3 (30) 



The energy spectrum is again anisotropic and may be estimated as, E(k±, k z ) k\ (A) 
w^, where k± ~ A -1 . Now, using equation (I3"0"j) . we have 



w 2 ^/3 



Summary of the balanced strong cascade 



1. The cascade is critically balanced in that the cascade time is of order of the wave period 
throughout the inertial-range. 
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2. The rms amplitudes have a scaling with the transverse scale which is of the Kolmogorov 
form. 

3. The turbulent cascade occurs in the transverse as well as parallel directions. However, the 
cascade occurs predominantly in the transverse directions, because the parallel correlation 
lengths scale as the 2/3rd power of the transverse correlation lengths. 



6. Imbalanced Weak MHD Turbulence 



The cascade is due to the resonant 3-wave interactions, see Galtier et al. (2000) and 
Lithwick & Goldreich (2003). Consider steady, imbalanced and isotropic excitation of ±ve 
Alfven waves, with rms amplitudes w~[ < <C Va and scale L, which initiates weak cas- 
cades to smaller spatial scales. As we have discussed earlier, the 3-wave resonance conditions 
forbid a parallel cascade, so wavepackets of transverse scale A < L have parallel scales L. 

The cascade times for the ±ve wavepackets can be estimated by considerations similar to 
those for the balanced case. Consider collisions suffered by a +ve wavepacket of transverse 
scale A and parallel scale L, with — ve wavepackets of similar spatial scales. A single colli- 
sion occurs over time ~ L/Va, during which the +ve wavepacket is perturbed by amount, 
Sw x ~ (Lw^w^ /XVa)- Successive perturbations add with random phases, so the number 
of collisions for perturbations to build up to order unity is N x ~ (XVa/ Lw^) 2 . Similar 
considerations apply to collisions suffered by a — ve wavepacket with many +ve wavepackets, 
and hence ~ (XVa/Lw^) 2 . Therefore the cascade times for the ±ve waves are, 



Kolmogorov's hypothesis of the A-independence of the energy fluxes 



K) : 



(33) 



imples that 



w^w x oc A (34) 

In other words, the requirement that both e + and e~ be independent of A leads to just one 
relation for the two quantities, w x and w^. 
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Hence scaling arguments are, by themselves, insufficient for the determination of the inertial- 
range spectra of the imbalanced cascade. The new physics that determines the spectra is the 
fact that the ± energies are forced to equalise at the dissipation scale; this pinning completely 
determines the spectra. 

Let the dissipation scale be at Adi S , and the common rms amplitude of and at this 
scale be u> dis . Then equation (1341) implies that we can write, 



wj = w dis 



^ \ (l+«)/2 



A 



dis 



A 



(l-a)/2 



W X = Wdis ( t — ) (35) 

Adis, 

where the parameters, Wdis, Adi S and a are to be determined. When the dissipation is caused 
by diffusive processes and v ~ r\ (i.e. the Prandtl number is close to unity), the dissipation 
timescale is 

tdis ~ ^ (36) 
u 

At the dissipation scale, the cascade time and the dissipation time are comparable. Setting 
A ~ Adi S in equation (1521) and setting it equal to £<jisj we can estimate that, 

T/ \ 1/2 



Note that w d i S is independent of the the energy fluxes, e ± . However, the other two param- 
eters, Adis and a depend on the energy fluxes, and these can be determined only by use of 
the kinetic equations of the weak turbulence theory of the cascade, as shown by Galtier et 
al. (2000) and Lithwick & Goldreich (2003)0. So, in contrast to the balanced weak cascade, 
dissipation plays a direct role in determing the inertial-range spectra of imbalanced weak 
MED turbulence. 

In common with the balanced weak case, we still need to verify that the turbulence is 
indeed weak in the inertial-range, So we need to verify that the quantitites, 



It should be noted that A^is and a can be expressed in terms of the rms amplitudes at the stirring scale, 
■f£ and w~^, but it is the energy fluxes, e^, that are the more interesting physical quantities. 
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Xi ^ « 1 (38) 

throughout the inert ial-range, L > A > A*. To do this, we use w^w^/X = u^ is /Adi S and 
equation ( 1371) to first work out Adi S : 

Adis ~ (39) 
wlw L 

For the inertial-range to be of large enough extent we must have L ^> Adi S , which implies 
that 



w T w T v , . 

Now, requiring that xf ^ 1 when A ~ A* gives us another condition: 



"' 7 "' 7 « (^r) 1 ' 2 (41) 



r2 



VI \LV A 

When v -C LVa, the above inequalities can be satisfied with appropriately chosen w^w^. 
and an imbalanced weak cascade with a substantial inertial-range is possible. 



7. Imbalanced Strong MHD Turbulence 

This is clearly the more general and generic case of MHD turbulence, and the account 
below is based on Lithwick, Goldreich & Sridhar (2007). Consider steady, imbalanced and 
isotropic excitation of ±ve waves, with rms amplitudes w~[ < < Va and scale L. 

At some transverse scale A C L, let the amplitudes be > w^, and parallel correlation 
lengths be A^ and A^". Note that we allow the parallel correlation lengths to be functions 
of the transverse scale, because resonance conditions are not in force in strong turbulence. 
+ve waves bend field lines by angle (w^/Va), hence a +ve wavepacket of parallel scale 
A^ will displace field lines in the transverse directions by distances {w^A^/Va)- If this 
distance exceeds the transverse correlation length, A, then the — ve waves will suffer strong 
perturbations. In other words, the — ve cascade is strong if 
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Then the — ve wave cascade time is 



+ AM > j (42) 
AA AV. - v ; 



t- ~ A (43) 



and the — ve wave parallel correlation length is 



A A ~ V A r x ~ A ( ^ ) (44) 



W A 



From equations ( 1421) and ( 1441) . we note that 



A+ > A- (45) 

We can infer a deeper result. Any two regions in the +ve packet that are separated by 
parallel lengths > A^ will be cascaded by statistically uncorrelated — ve waves. Hence the 
— ve waves will imprint their parallel scales on the +ve waves, and both ±ve waves will have 
the same parallel correlation lengths, A\ : 

A+ ~ A- ~ A A ~ A (J^J (46) 



From equations (1421) and (1461) we have 



Xl ~ 1 (47) 



Therefore the —ve cascade is critically balanced. 
It turns out that the +ve wave cascade time is, 



^A 



This means that the straining rate imposed by the — ve waves on the +ve waves, (w x /A), is 
imposed coherently over time (X/w x ). Then 
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— ve wave coherence time ttVA wt , , 

> 1 4 9 

— ve wave wave period Aa w x 

How can the coherence time of the — ve waves exceed their wave period? The answer is 
involved and is worked out in the Appendix of Lithwick, Goldreich & Sridhar (2007); here 
we merely summarise the result. Now let us look at the interaction between ±ve waves of 
transverse scale A and parallel scale Aa, in the rest frame of the +ve waves. In this frame, the 
— ve waves alter +ve waves which react back onto the — ve waves. Let z' measure distance 
along the parallel direction in the rest frame of the +ve waves. At a fixed z' location the 
+ve waves are changing on their cascade time scale t~^. Hence, over times separated by t^, 
the — ve waves crossing z' are cascaded by entirely different +ve waves. This implies that 
t~ rrA ~ t^. Because the +ve waves are strained at rate (w^/X), it follows that t x ~ (X/w^). 

We are now in a position to work out the scalings for the imbalanced strong cascade. 
Kolmogorov's hypothesis of the A-independence of the energy fluxes 



A 



+ \2 

(50) 



t\ X 



implies that in the inertial range, 



± 



A 1 / 3 



(e=F)i/3 



A ^ ~ ^ K,A (51) 

Equations (150]) are given in Verma et al. (1996) and Verma (2004), although isotropy is 
assumed. 

Summary of the imbalanced strong cascade 

1. The ratio of the Elsasser amplitudes is independent of scale, and is equal to the ratio of 
the corresponding energy fluxes. Thus we can infer turbulent flux ratios from the amplitude 
ratios, thus providing insight into the origin of the turbulence. 
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2. In common with the balanced strong cascade, the energy spectra of both Elsasser waves 
are of the anisotropic Kolmogorov form, with their parallel correlation lengths equal to each 
other on all scales, and proportional to the 2/3rd power of the transverse correlation length. 

3. The equality of cascade time and waveperiod (critical balance) that characterizes the 
strong balanced cascade does not apply to the Elsasser field with the larger amplitude. 
Instead, the more general criterion that always applies to both Elsasser fields is that the 
cascade time is equal to the correlation time of the straining imposed by oppositely-directed 
waves. 

4. In the limit that the energy fluxes are equal, the turbulence corresponds to the balanced 
strong cascade. 

We note that there are other theories of imbalanced cascades which differ from the 
above account; to the best of my knowledge, these are given in Beresnyak & Lazarian (2008), 
Chandran (2008), Perez & Boldyrev (2009) and Podesta & Bhattacharjee (2009). 



8. Stirring-affected scales in the imbalanced strong cascade 

Consider again steady, imbalanced and isotropic excitation of ±ve waves, with rms 
amplitudes < < Va, on scale L. The wave periods of the ±ve waves at the stirring 
scale L/Va- Let the coherence time of the waves on the stirring scale be T; we expect 
that T > L/Va, but otherwise T is unconstrained and determined only by the physics of 
however the ±ve waves are generated on scales ~ L. The theory of Lithwick, Goldreich & 
Sridhar (2007) applies when T > t\\ otherwise, when T <t~£, the coherent straining of +ve 
waves by the — ve waves will be interrupted before the +ve waves can cascade. Using, 

*t ~ ^corrA ~ ^ " (52) 

we can infer that T <t^ implies, 

{£-) T z l 2 

1 = ~^yJT < L (53) 

Then, the results of Lithwick, Goldreich & Sridhar (2007) are valid for only for transverse 
scales A < I. 

The range of transverse scales, £ < A < L, can be referred to as stirring-affected scales. 
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For these scales, the +ve wave cascade is no longer strong and ±ve spectra are, as yet, 
unknown. It would also be interesting to see if there is any signature of these stirring- 
affected scales in solar wind turbulence. 



9. Some comments 

We have assumed that the strong large-scale magnetic field is uniform in magnitude and 
direction, but this need not be the Kraichnan (1965) realised. We could consider a 

more general setting in which the magnetic field is disordered on large scales, with the largest 
spatial correlation length, L, being smaller than the system size. Let the magnetic energy 
in eddies of size L be larger than the magnetic energy on smaller scales. Then the magnetic 
field of the eddies of size L will act as the mean magnetic field for Alfven and Slow waves of 
smaller scales, and all of what we have derived in this review for the inertial-range of MHD 
turbulence in a strong magnetic field may be expected to be still validE The opposite case 
of the large-scale magnetic field being weak often occurs in problems concerning dynamo 
processes, and turns out to be harder to analyse. 

We have considered incompressible MHD turbulence, but the general case in astro- 
physics involves compressive fluids. In compressible MHD, there are four types of waves: 
Alfven, Fast, Slow and Entropy. The inertial-range spectra we have derived are applicable 
to cascades of Alfven waves, which are always incompressible. These are valid, so long as the 
other three waves (which are compressive) have negligible effect on the Alfven wave. That 
this is indeed so has been discussed in the following papers: Lithwick & Goldreich (2001), 
Maron & Goldreich (2001) and Cho & Lazarian (2003). 

In the context of Navier-Stokes turbulence of a neutral fluid, it is well-known that a 
triple correlation function of velocities is related to the rate of energy input; see Karman & 
Howarth (1938), Kolmogorov (1941) and Monin (1959). Similar relations in the context of 
MHD turbulence have been obtained by Politano & Pouquet (1998a,b) and Podesta (2008). 
However, the relation to the theories of the inertial-range spectra of MHD turbulence has 



5 There is an interesting detail here that is, perhaps, worth noting. In the case of hydrodynamic turbulence, 
the velocity field of the large-scale eddies (of size L) will "sweep" eddies of smaller sizes, and this effect gives 
rise to observational consequences for the velocity time correlation function measured by a probe fixed at 
some location in the fluid. However, sweeping cannot affect the dynamics of nonlinear interactions, because 
a large-scale velocity field can be effectively transformed away by an appropriate Galilean boost. On the 
other hand, a large-scale magnetic field cannot be similarly transformed away; one way to physically see this 
is to note that such a field enables (Alfven and Slow) waves that travel both up and down the field lines, 
and these cannot be wished away by any Galilean boost. 
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not received as much attention as in the case of neutral fluids. 

It has sometimes been argued that there is a "dynamical alignment" of velocity and 
magnetic fields, resulting in spectra that are flatter than Kolmogorov — see Boldyrev (2005) 
and Beresnyak & Lazarian (2006). However, observations of solar wind turbulence show 
spectra that are consistent with Kolmogorov - see Horbury, Forman & Oughton (2005). 

I would like to thank (i) the organisers of the meeting on Astrophysical Magnetohydro- 
dynamics for their hospitality, (ii) Peter Goldreich and Yoram Lithwick for all the things I 
learnt from them. 
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